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Polymer quantization is as a useful toy model for the mathematical aspects of loop quantum 
gravity and is interesting in its own right. Analyzing entropies in the standard Hilbert space and 
the polymer Hilbert space we show that they converge in the limit of vanishing polymer scale. We 
derive a general bound that relates entropies of physically equivalent states in unitarily inequivalent 
representations. 
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Statistical thermodynamics, quantum mechanics and 
information theory promoted entropy from an auxiliary 
variable of the mechanical theory of heat [l| to one of the 
most important quantities in science 0, [J. Microscopic 
derivations of the thermodynamical entropy were, and 
still are, a major part of the statistical physics. Entropy 
and related quantities underpin most of the classical and 
quantum information theory Black hole thermody- 
namics Q introduced the idea that an upper bound on 
the entropy of a system (and entanglement between the 
subsystems) scales with the area, and not with the vol- 
ume. Formalized as a holographic principle, the area law 
was demonstrated in a variety of situations @, @] , which 
include string and loop quantum gravity (LQG) Q cal- 
culations of black hole entropy @,@] and related models 

LQG and the spin networks/ spin foams formalisms 
represent the idea of fundamental discreteness. They 
are conceptually and technically different from the usual 
quantum field theory on a space-time continuum. This 
difference leads to a number of questions. One of them 
@) HH is emergence of the semiclassical states of quan- 
tum gravity. These states have to enable the space-time 
continuum view and the ensuing construction of a field 
theory. Investigation of this problem and the associated 
mathematical issues were one of the motivation to study 
the polymer quantization of a single non-relativistic par- 
ticle [ll| as a useful toy model. 

Another question pertains to the validity of holo- 
graphic principle(s) and compatibility of the entropy cal- 
culations across different quantization schemes. Recent 
advances in coupling gravity and matter [l2| in LQG [1] 
make it necessary to re-consider holography in these sys- 
tems, which is derived separately in the gravity and the 
matter sectors. One can wonder if two unitarily inequiva- 
lent quantization schemes give compatible predictions for 
observables, are the predicted entropies "close"? While 
the desirable answer is "yes", its validity is neither ob- 
vious nor guaranteed. Again, the polymer quantization 
provides a convenient model setting to study this prob- 
lem. 

We demonstrate that, if the Schrodinger and the poly- 
mer schemes give "close" (to be precisely defined be- 
low) predictions for fundamental observables, their pre- 
dictions of entropy are close. The two entropies coincide 
in the continuum limit of the polymer quantization. We 
illustrate this by comparing the entanglement entropy in 



the ground state of coupled harmonic oscillators in the 
two schemes. 

The rest of this paper is organized as follows. After 
a brief review of the polymer quantization we describe 
our main result, outline the calculational techniques for 
Gaussian states, and present the example. We conclude 
with discussing possible limitations and future directions. 



I. REVIEW OF POLYMER QUANTIZATION 

A useful starting point is a unitary Weyl (displace- 
ment) operator, 



W(q,k)^ = W(-q,-k), 



(1) 



where q and k are the translation an boost parameters, 
respectively. The Weyl operator can be expressed as a 
product of the translation and boost operators, 

W(q, k) = ei qk V{q)U{k) = W(q, k) = e -? qk U(k)V(q), 

(2) 

Here and in the following we set h = 1 . 

In the standard Schrodinger representation these uni- 
tary operators are generated by the self-adjoint position 
and momentum operators, 



V(q) = e lqp , U(k) = e lkx , (3) 

and the Weyl operator acts on wave functions as 

W(q, k)tp(x) = e~^ kq e ikx ^{x + q). (4) 

The polymer representation is built around (Harald) 
Bohr compactification of the real line [l4|, EH • A slightly 
different construction is possible by noting that the (po- 
sition) space Hpoiy is spanned by the basis states \x), 



ipx(x') = (x'\x) = (x'\x) = 5 X 



(5) 



where S is Kronickcr's delta, and the various bracket signs 
are explained below. Countable superpositions of these 
states, 



I*} 



ip(x)\x), 



(6) 



where 7 is a countable set (a "graph" [llj ) and functions 
-0 satisfy certain fall-off conditions, form a space of cylin- 
drical functions Cyl 7 . The infinite-dimensional space Cyl 
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is a space of functions that are cylindrical with respect to 
some graph. The Hermitian inner product on Cyl follows 
from Eq. ©, 



(7) 



The Cauchy completion of Cyl is Hpoiy, and the the 
triplet of spaces 



Cyl C H po i y C Cyl*, 



(8) 



where Cyl* is an (algebraic) dual of Cyl, shares simi- 
larities with the construction of the physical space of 
LOG llf . The questions of topology are discussed in 
[ill llZL 03 ■ Dual elements are denoted by (T|, and their 
action on the elements of "Hpoiy by (T^) = T(^). The 
inner product defines a dual element by 

$(#) = ($|#). (9) 

We introduce the momentum dual states (p\ £ Cyl* by 

(p\x) = e~ ipx . (10) 

Hence 

(p\ = j2 e ~ ipx ( x i m = (pm- (11) 

In the "position representation" that was described so 
far the operator U(k) has a self-adjoint generator x, that 
acts by multiplication and has normalized eigenstates, 

x\x)=x\x). (12) 

However, since \xj) and V(q)\xj) are orthogonal to each 
other no matter how small the translation is, the transla- 
tion operator is not weakly continuous and the momen- 
tum operator is undefined. Similarly, while it is easy to 
show that the usual relation 



(p\V(q) = e**(p| 



(13) 



holds, it docs not correspond to any state of Cyl. 

On the other hand, a momentum representation is 
spanned by the states \kj). In this case there exists a self- 
adjoint momentum operator with normalizablc cigncs- 
tates, but the position operator is undefined. Since 
the polymer representation(s) are not weakly continuous, 
they do not satisfy the conditions of the Stone- von Neu- 
mann [T3I ] theorem and thus are not unitarily equivalent 
to the Schrodinger representation [III [l4| ■ 

The two constructions of the "H po iy allows the most 
extreme expression (l4j of Bohr's complementarity fltij |. 
Failure of the weak continuity makes it into a convenient 
toy model of LQG, and it was recently studied both as 
such and in its own right (ill [l7l - [l9j . Moreover, the poly- 
mer quantization and its generalization to a scalar field 
were applied to the problems of quantum cosmology and 
spherically-symmetric collapse [20(. 



Fell's theorem [2lJ guaranties that, although two rep- 
resentations may be unitarily inequivalent, by using a 
finite number of finite-precision expectation values of ob- 
servables it is impossible to distinguish between the two. 
More precisely, for a state p\ of one of the representa- 
tions, a finite set of operators Ai (belonging to the Weyl 
algebra) whose expectations are calculated on that state, 
and the set of tolerances there exists a "physically 
equivalent" state P2 on another representation resulting 
in the expectation values that differ from the first set 
by less than the prescribed tolerances. Since its proof is 
not constructive, an explicit construction of states and 
operators is needed in each case. 

While a mathematically rigorous construction involves 
a multi-scale lattice [13], we follow the approach of [TTJ] 
and will work with a regular lattice 7, where the neigh- 
boring points have a fixed spacing p. An effective mo- 
mentum operator is introduced through the finite differ- 
ence 

P^-^(V{p)-VHp)). (14) 

The limit p — > in the Schrodinger representation gives 
the usual momentum operator p. Similarly, its square is 

Pi = (Pj 2 = i ( 2 - ^) - ^(-2/0) (15) 

The operators are self-adjoint fli1.li~7l]. However, p^ is not 
positive on a general Cy L e ven if it is a regular lattice. 
For example, for a state [l7| 



(16) 



££7 



the expectation value (^Ip^l^) is real only for a sym- 
metric graph, i.e. x £ 7 —x £ 7. We will restrict 
ourselves only to such graphs. 

It follows from Eq (jT5j) that p^ skips the neighboring 
lattice sights when acting on states, hence its eignefunc- 
tions can have support on either even- or odd-numbered 
sights n/i, n £ Z. As a result, using p 2 ^ as the kinetic 

terms for a Hamiltonian H a leads to a double degener- 
acy of the eigenstates [III [TJl when compared with the 
Schrodinger representation. Without appealing to it, one 
can note that the state counting gives twice the semi- 
classical result fdpdx/2ir. Depending on one's goals it 
is possible either to adjust the state counting or to intro- 
duce a kinetic operator [llj] 

K = P\ = ( 2 - - V(-p)) , (17) 
The commutation relations between the operators are 



[x 1 V( P )]=-pV(ii), 



(18) 



[*>P»] = o(^G«) + v(-n)) = i(i - ±p 2 kj 7 (19) 



(20) 
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A key step in extracting physical predictions in the 
polymer quantization is to consider the shadow states 
fill ], that realize the physically equivalent state on W po i y . 
These are projections of the elements of Cyl* onto Cyl 7 , 

($i = ^2§xx)(x\ -> i$ shad ) = J2 ®( x )\ x )- ( 2i ) 



xEM. 



xg7 



Using these states it is possible to explicitly demon- 
strate physical indistinguishability of the predictions of 
the polymer and Schrodinger quantizations, as mandated 
by Fell's theorem. 



II. ENTROPY 

Unlike observables that can be constructed from the 
Weyl algebra, there is no existence result that guaranties 
equivalence of the predictions for entropy. Moreover, it is 
known 0] that entropy is not a continuous function, and 
without additional restrictions there are states of infinite 
entropy in a neighborhood of any state. 

Consider a Schrodinger state 



(22) 



which can be built from the states j) by the Gramm- 
Schmidt procedure. As the lattice scale p decreases, the 
shadow states tend to an orthonormal system, 



(26) 



hence if the basis vector 6 /C M is obtained from 
|*(/i)i>, then 

l*0*)k> =£<»>, I ™ C fc l =*«- ( 27 ) 

All separable Hilbcrt spaces are isometrically isomor- 
phic, so for every p we can introduce a unitary operator 
Ufj, : ICfj, — > /C, where K, is some "standard" Hilbert space. 
Then 

U»%) = \1)ZK, Cg*(M) fc >=5>£|Z>. (28) 



Consider two density operators on the space IC, 

p = J2^i\i)(i\, ^ = £™i«IW (29) 

i ikl 

where S(p) = S(p) and S(p po i y ) = S{p fl ). Using Eq. (HTJ) 
we find that 



where the form an orthogonal basis. If this basis is 
formed by the eigenstates of some operator, and its poly- 
mer counterpart has the corresponding eigenbasis 
then the decomposition 



n° 

Fpoly 



(23) 



trivially has the same entropy, but it is not guaranteed 
that the expectation variables of all the observables of 
interest are close. 
The shadow analog 



Ppoly 



(24) 



with \ty(p)i) being the normalized shadow states 
|^ r (/i)| had ), that are close "approximations" of the states 
in the sense of Fell's theorem. The shadow states 
are not in general orthogonal [11], making the entropy 
estimation more complicated. 

Concavity of entropy Q leads to an inequality that in 
our case reads as 



SVpoly) < - Wj log i 



S( P ), 



(25) 



This result holds for any two representations were the 
equivalent state of any is pure. 

Consider a (completion of the) linear span of the 
shadow states \$>(p)i) as a Hilbert space JC^. It is a 
separable Hilbert space and has an orthonormal basis, 



lim tr|,o - pj = 0. 



(30) 



The lower semicontinuity of entropy Q gives 

S(p) <]iminf£(p M ), (31) 

which together with the bound (|25[) establishes the con- 
vergence 



lim SQspoiy) = S(p). 



(32) 



III. COUPLED HARMONIC OSCILLATORS 

Wc illustrate our result by considering entanglement of 
the ground state of two oscillators with the Hamiltonian 

H = ^(Pl+Pi) + 4) + A(*i - x 2 )\ (33) 

A parameter that determines the closeness of the physical 
prediction is the ratio of the lattice size p to the oscillator 
scale 



[moj)- 1 ' 2 . 



(34) 



The ground state of this Hamiltonian (in the 
Schrodinger quantization) is Gaussian [2^, HU, i.e., all 
statistical moments of the canonical observables can be 
expressed from the first two moments — the expectation 
values and the variances. Ground states of any quadratic 
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n-particlc Hamiltonian, as well as thermal states, coher- 
ent states and squeezed states are of this type. Gaussian 
states are very important in quantum optics and quan- 
tum information processing. They have a number of use- 
ful mathematical properties, of which we need two. First, 
a reduced density matrix of a Gaussian state (e.g., den- 
sity operator of the first oscillator in the above example) 
is also Gaussian. Second, entropy of a Gaussian state can 
be calculated using its symplectic eigenvalues, as we now 
describe @,[H. 

Position and momenta (either classical or quantum) of 
n particles can be arranged into a single 2n-dimensional 
vector r T = (xi, . . . , x n , p\, . . . , p n ). Classical Poisson 
brackets and quantum commutation relations can be ex- 
pressed using the 2n x 2n symplectic matrix J, 



m 



(40) 

The correlation term vanishes exactly. Hence, keeping 
only the leading terms in p/d the correlation matrix be- 
comes 



(d 2 o 



(41) 



I n 
-In 



The product of uncertainities the pol ym er quantisation is 
less than its standard value of j [n], [24[ . The state is no 
longer exactly Gaussian: because the correlation matrix 
violates the inequality T u + i J > 0, and Eq. ([55)) gives a 



where /„ is the n-dimensional identity matrix, and the 
symplectic matrix satisfies J T = —J = J -1 . The Hamil- 
ton equations are given by r = JdH/dr, and the canon- 
ical transformations are represented by matrices S that 
satisfy SJS T = J, thus forming the n-dimensional sym- 
plectic group. 

The vector of expectation values and the symmetric 
correlation matrix are defined by 



[Vfc,rfc] = iJu, (35) complex value of the entropy, S ~ £iL. s for a pure state. 

The standard measure of a pure state bipartite en- 
tanglement is the von Neumann entropy of either of the 
reduced density matrices Q. For the Gaussian states 
it can be calculated using Eq. (|38|) . whith symplectic 
eigenvalues for the correlation matrix of the subsystem. 
Tranforming to the normal modes (this is a simplectic 
transformation) gives two uncoupled oscillators with the 
frequencies 



D 



{Wi ~ A), (fj - Dj)}) (36) 



where {a, 6} is the anticommutator. Unitary operators 
U(S) that by transforming the Gaussian state according 
to p' = U(S)pW (S) transform the statistical moments 
according to 



D' = SD, 



V = STS 1 



(37) 



are called Gaussian operators. They represent important 
experimental procedures in quantum optics and quantum 
information. 

Any Gaussian correlation matrix T can be di- 
agonalized by some symplectic transformation Syy, 
Sw^S^ = W, where W is a diagonal matrix 
diag(ci, . . . , <7„, a%, . . . , a n ). This spectrum consists of 
the absolute values of the eigenvalues of the matrix i JT. 
The eigenvalues crj are the symplectic eigenvalues of the 
covariance matrix. Finally, entropy of a n-particle Gaus- 
sian state p is given by its symplectic eigenvalues as [7L [23j] 



UJ 2 = OJ2 = W\/ 1 
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(42) 



The (Schrodinger) correlation matrix in the normal co- 
ordinates is r' = diag(c? 2 , d 2 /a, d~ 2 , d~ 2 a), and the only 
eigenvalue of the reduced correlation matrix is 



ci = 



1 + a 
V5" 



(43) 



The polymer quantization is treated similarly. Assum- 
ing a > 1 and keeping only the leading order terms in 
the powers of p/d we find 



T p°iy 



(1 + a 2 ) ^ 2 

d 2 



01 



(44) 



and 



S^ = S- l -±flo g ^^. (45) 
y/a — 1 d A 



3=1 V 

(38) 

Covariance matrix of the ground state of a single har- 
monic oscillator (in the Schrodinger representation) is 

/ 2(£ 2 ) (xp) + (px)\ _(d 2 0\ , , 
1 - ^ (px) + (xp) 2{p 2 ) J ~\Q d- 2 ) ■ [ ™> 

Calculations in the polymer representation also give 
zero expectations, (x) = 0, (p n ) = 0, but the variances 



IV. SUMMARY AND OUTLOOK 

We showed that the shadow state construction gives 
not only "close" expectation values for the observablcs 
in the unitarily inequivalent representations, but the en- 
tropies of the two states agree as well. The lattice ef- 
fects modify the expectation values and commutation re- 
lations. Shadow states do not satisfy the Gaussian prop- 
erty exactly, but only up to the terms of the order of 
p 2 /d 2 . 
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Since the convergence of entropy of the shadow states 
to the Schrodinger representation value is not necessarily 
uniform, the following scenario is plausible. Both the dis- 
crepancy in the expectation of the momentum variances 
and symplectic eigenvalues are of the order of n 2 /d 2 for 
each (uncoupled) oscillator. This is also order of magni- 
tude of the corresponding change in entropy contribution 
if Eq. is used. Hence, even if the expectations of 
the observables agree, for a fixed value of \x and a large 
number n of oscillators the two entropies will differ by 
~ nji 2 /d 2 , which may be a significant amount. We leave 



for future work a precise estimate of the discrepancy in 
entropies for a fixed scale and a large number of particles. 
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